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ABSTRACT 


Context. Supermassive stars and quasi-stars (massive stars with a central black hole) are both considered as potential 
progenitors for the formation of supermassive black holes. They are expected to form from rapidly accreting protostars 
in massive primordial halos. 

Aims. We explore how long rapidly accreting protostars remain on the Hayashi track, implying large protostellar radii 
and weak accretion luminosity feedback. We assess the potential role of energy production in the nuclear core, and 
determine what regulates the evolution of such protostars into quasi-stars or supermassive stars. 

Methods. We follow the contraction of characteristic mass scales in rapidly accreting protostars, and infer the timescales 
for them to reach nuclear densities. We compare the characteristic timescales for nuclear burning with those for which 
the extended protostellar envelope can be maintained. 

Results. We find that the extended envelope can be maintained up to protostellar masses of 3.6 x 10° m3 Mo, where 
m denotes the accretion rate in solar masses per year. We expect the nuclear core to exhaust its hydrogen content 
in 7 x 10° yr. If accretion rates rh >> 0.14 can still be maintained at this point, a black hole may form within the 
accreting envelope, leading to a quasi-star. Alternatively, the accreting object will gravitationally contract to become 
a main-sequence supermassive star. 

Conclusions. Due to the limited gas reservoir in typical 10 Mo dark matter halos, the accretion rate onto the central 
object may drop at late times, implying the formation of supermassive stars as the typical outcome of direct collapse. 
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. Introduction 


upermassive black holes (SMBHs) with more than 10° Mo 
have been observed at z > 6 (Fan et al! 
| 2006), and recently even at z 7.085 
(Mortlock et al. . While theoretical considerations 
indicate a variety of potential pathways, including direct 
collapse to a black hole, collapse to a supermassive star 
collapse of stellar clusters or clusters of black holes 
[Volonteri & Bellovary|[2012), substan- 
tial efforts have been made to assess the validity of such sce- 
narios. As the most straightforward possibility, one could 
imagine these SMBHs to originate from remnants of the 


first stars, which were proposed to be very massive with 
100 — 1000 Mo (Abel et all 0 Brom Taron 2004). 
Stellar evolution calculations by (2002 
indeed show that black holes may form for stellar masses 
between 10 — 30 Mo, or 100 — 300 Mo, thus yielding the 
seeds for further growth to SMBHs. 


New calculations following the evolution of primor- 
dial star formation beyond the first peak however show 


that fragmentation can occur efficiently (Stacy et al. 
; ; (Greif et_alJ|2011; 


the accretion luminosity on the protostellar core seems to 


However, if high accretion rates are maintained, a quasi-star with an interior black hole may form. 


be minor (Smith et_al.|2011|, , UV feedback during the 


main sequence phase appears to set the upper limit to the 


stellar mass to 50 — 100 Mo (Hosokawa et al. Bon, (20120; 
[Susal2013). Even if very massive stars could form, the re- 
sulting HII region is likely to inhibit accretion for at least 
108 yrs Moai ié ot apod). An additional prob- 
lem has been identified by (2012), who 
showed that black holes resulting from such stars would 
be born with strong dynamical kicks, implying that they 
should readily be expelled from the halo. 

On the other hand, dense stellar clusters may form 


at high redshift in the presence of trace amounts 
of dust (Schneider et al) [2003 2005 


(Clark et _al]|2008; [Omukai et_al] [2008 2011; 
Schneider et al] |2011; [Dopcke et al! 
2012; 2012; . The relativistic 
instability of such clusters_has been explored in a semi- 
analytic framework by and 
(2010, , finding characteristic black 
hole masses of up to ~ 3000 Mg. The latter is already sub- 
stantially higher than the masses of the first stars, though 
it may still be difficult to obtain masses as high as 10? Mo 
at z = 7. In the presence of a significant black hole spin, 
which will be rapidly obtained in the presence of efficient 
accretion, one typically expects an additional growth by 
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about four orders of magnitude, thus requiring seed masses 
of the order 10° Mo {Shapird 2005). 

The formation of such massive seeds has been con- 
sidered in the context of the direct collapse model, 


where gas in a 107-108 Mo halo is expected to 
collapse without fragmentation _into ingle 


order to avoid fragmentation, molecular hydrogen needs 
to be efficiently dissociated, requiring a strong am- 
: 2010 2011; 

). Such radiation backgrounds 
may have been provided by nearby starburst galaxies, which 
may indeed occur frequently enough to explain the ob- 
served abundance of SMBHs at z ~ 6 (Didkstra et al/fo008 
Agarwal et_al][2012). 

Using numerical simulations, have 
modeled the gravitational collapse of massive primordial 
halos cooling via atomic hydrogen, reporting an isother- 
mal density profile and angular momentum transport via 
bar-like instabilities during the formation of the first peak. 
presented the first study follow- 
ing the evolution beyond the first peak, and modeling the 
formation of self-gravitating disks on parsec-scales. While 
these studies employed a typical resolution of 16 cells per 
Jeans length, it was recently demonstrated that turbulent 
structures can only be resolved with a resolution of at least 
32 cells per Jeans length, preferably more 
. Ina set of high- 
resolution simulations following the formation of the first 
peak, [Latif et al] (2013allh) found extended turbulent struc- 
tures in the center of massive primordial halos, but no 
signs of simple disks or bar-like instabilities. Disks do how- 
ever form at later stages of the evolution, with masses of 
~ 1000 Mo on scales of 30 AU, and characteristic accretion 
rates of ~ 1 Mg yr7 (Latif et alJ2013¢). As a result, very 
massive central objects may indeed form in rather short 
cosmic times. The dynamics at early times have also been 
explored by (2013), finding a somewhat differ- 
ent result including the formation of toroidal structures on 
scales of several parsec. It is however not fully clear to 
which extent the latter is a result of the simplified initial 
conditions employed in their calculation. In order to pro- 
vide more quantitative predictions regarding the conditions 
where black holes may form, explored 
the correlations of the baryon spin with the spin of the 
dark matter, showing that the resulting correlation is how- 
ever weak, and that the baryon properties cannot be naively 
extrapolated from the dark matter. 

While the accumulation of high masses seems feasible 
from a hydrodynamical point of view, the resulting object 
is still considerably more uncertain. As a first step, we 
expect the formation of a massive protostars, as the gas 
becomes optically thick at high densities 
. It is however unclear 
how these objects are going to evolve, and whether they 
will form a supermassive star or a quasi-star. While a 
supermassive star denotes a conventional star with very 


high masses of 10° —10° Mo (Shapiro & Teukolskyi1986), a 


quasi-star refers to an object with similar mass, but where 
the central core has collapsed into a low-mass black hole 
(Begelman et_al.| |2006: [Begelman 2010). The latter then 


Article number, page 2 of [9] 


accretes mass from the stellar envelope, while the quasi-star 
as a whole may accrete at a rate larger than the Eddington 
rate of the black hole. A central question however concerns 
the conditions under which we expect the formation of a 
quasi-star as opposed to a supermassive star, and which of 
these objects should be expected as the generic outcome in 
the context of the direct collapse model. Also the proper- 
ties of these objects are of high interest, as the amount of 
stellar feedback may influence the final accretion rates. 

In_ order to explore that question, 
recently followed the evolution of rapidly accret- 
ing protostars, showing that they expand as cool super- 
giants, thus inhibiting the feedback from accretion lumi- 
nosity and potentially allowing accretion to proceed for 
very long times. Hosokawa et al. also found that nuclear 
burning starts at stellar masses of about ~ 50 Mo, and is 
maintained until the end of their calculation at 10° Mo. 
Therefore, their model indicates no transition towards a 
quasi-star at least during the evolutionary phase consid- 
ered. It is however important to assess how long this phase 
of efficient accretion can be maintained, and under which 
conditions a transition to a supermassive main-sequence 
star can be expected, which can potentially influence the 
accretion flow via radiative feedback 
[2002s (2012). While the majority of 
such supermassive stars directly collapses into a black hole 
(Fryer & Heger Bori), a small mass window exists around 
55000 Me where violent supernova explosions of up to 
10°° erg may occur (Johnson et al.|(2013). 

Following a different approach, 
and have proposed the formation of 
a black hole in the interior of rapidly accreting ob- 
jects, leading to quasi-stars as the progenitors of SMBHs. 
(2006) argued that for quasi-stars with 
10*-10° Mo, the typical accretion timescale is considerably 
shorter than the timescale for nuclear burning, and as a re- 
sult, the latter may not be able to stop the collapse of the 
central core. A potential advantage of such a configuration 
is that the accretion onto the central object is not limited 
by the Eddington accretion rate of the black hole, but by 
the Eddington accretion rate of the more massive quasi-star 
(Begelman et al/2008). considers objects 
with more than 10° Mo, which first evolve towards the 
main sequence as a supermassive star, but then form a black 
hole after an extended phase of nuclear burning. |Ball et_al. 
followed the evolution of black holes in a quasi-star 
employing the Cambridge STARS stellar evolution pack- 
a [Pols et alll190), finding that the black 
hole may efficiently accrete 10% of the stellar mass before 
hydrostatic equilibrium breaks down. They also report that 
the results are sensitive to the boundary condition in the in- 
terior. The contraction of an embedded isothermal core in 
a stellar envelope was further explored by [Ball et al] 
in the framework of the Schönberg-Chandrasekhar limit. 

While the models for the quasi-stars typically employed 
a Thompson-scattering opacity, re- 
ported that H~ rather than Thomson scattering dominates 
the opacity in the protostellar atmosphere. They find lu- 
minosities close to the Eddington luminosity, 


4nGM myc 4 M 
— =3. 1075 — |, 1 
a 3.8 x 10 (Z£) (1) 


where G is the gravitational constant, mp the proton mass, 
c the speed of light, or the Thomson scattering cross sec- 


Lgaa = 
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tion and M the mass of the star. The temperatures in the 

atmosphere are rather cool, ~ 5000 K, as the protostars are 

on the Hayashi track. As a result, the characteristic radii 
evolve as 

1/2 

M 

Rini = 2.6 x 10° Ro (=) . (2) 


The protostars are thus considerably more extended and 
show an increasing stellar radius as a function _of mass 
while the models of [Begelman et al] and 
indicate a constant radius as a function of stellar 
mass. The difference is crucial, as the stellar radius regu- 
lates the temperature on the surface, and thus the strength 
of protostellar feedback. In fact, the behavior reported by 
is well-known in cases where the 


accretion timescale, 


(3) 
is much smaller than the Kelvin-Helmholtz timescale, 


GM? 


However, report that this behav- 
ior extends into the regime where txy < tacc, i.e. where 
Kelvin-Helmholtz contraction is faster than mass growth 
via accretion. The physical mechanism which allows this 
phase to continue will be discussed in detail in section 2] 
Taking the maximum luminosity of the star 


M 11/2 R =1/2 
tmae~08b0 (a5) (Ge) k 


which assumes Kramer’s opacity k œx  pT~3/5 


(Hayashi et al| [1962) and an initial radius as given in 
Eq. (2), they show that both timescales are equal at 
protostellar masses of 


tKH = 


a 0.26 
Meg = 14.9 Mo (=) (6) 


where we parametrized the accretion rate as M = 
rn Mo yr t. On the other hand, the accreting star remains 
on the Hayashi track up to stellar masses of 1000 Mo, the 
highest mass reached in the calculation of 
20123). The extended envelope locked at Teg ~ 5000 K 
allows efficient accretion with only moderate feedback for a 
longer period. 

In this paper, we aim to assess how long this efficient 
accretion phase can be maintained beyond M ~ 1000 Me 
without strong feedback from the protostar. In addition, we 
will discuss the potential impact of nuclear burning for the 
evolution of the accreting objects, and the conditions under 
which a quasi-star as opposed to a supermassive star may 
form. For this purpose, we explore the interplay of mass 
accretion and Kelvin-Helmholtz contraction in section 
We calculate the impact of nuclear burning in section [B] 
A final discussion of our results is provided in section [4] 
The impact of additional processes such as deuterium shell 
burning and hydrogen burning via the pp-chain is assessed 
in the appendix. 
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Fig. 1. The enclosed mass as a function of radius in the 
protostar calculated from Eq. (10), assuming a total stellar mass 
of 10° Mo and mass accretion rates rn = 0.1, 1., 10. The current 
time t is calculated as t = (M/Mo)/m. For a given mass shell, 
larger accretion rates imply larger radii, as the star had less time 
to contract. 


2. The interplay of mass accretion and 
Kelvin-Helmholtz contraction 


In the following, we will sketch why the protostar may 
maintain large envelopes even if the Kelvin-Helmholtz 
timescale becomes shorter than the timescale for accretion. 
For this purpose, we consider mass shells of enclosed mass 
M, located at radii R(M,t). A mass shell M forms at the 
time M/M, with an initial radius as given by Eq. (2). The 
radii of these mass shells evolve on their Kelvin-Helmholtz 
timescale, given as 


GM? 


ne M D= 
xu(M, R) R(M, t)Luaa 


(7) 


As a result, we have 


dR R(M,t) 


dt ~~ tgu(M,R) (8) 


The equation is integrated between the initial radius Rini 
when the mass shell forms, given in (2), corresponding to 
the time tinj(M) = M/M, to the current radius R at time 
t. The integration yields 


1 1 
R 7 Rini 


ATMyC 
op M 


(t — tins(M)). (9) 


Inserting the initial radius as given in Eq. (2) and expressing 
the result in astrophysical units, we obtain 


1000Ro 1000 Ro 


1.04 yr (t = tini(M)) 
R 2.6 x 102(M/Mo)1/2 


M/Mo 

(10) 
The first term on the right-hand side thus dominates right 
after the formation of a mass shell and determines its initial 


radius, while the second term describes its evolution due to 
Kelvin-Helmholtz contraction. It is remarkable that at late 
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times, the second term dominates and the evolution of the 
radii appears to be independent of their initial position. 
This can be understood, as our expression for the Kelvin- 
Helmholtz timescale, Eq. (Z), scales with the inverse radius 
of the mass shell. As a result, the evolution slows down 
during the contraction, implying that the Kelvin-Helmholtz 
timescale may increase significantly in the interior. An ex- 
ample for the resulting structure in a 10° Mo star is given 
in Fig. [] including both the regime where the first and the 
second term of the equation dominate. 

We now aim to estimate when a given mass shell reaches 
the densities of nuclear burning, which we take as pnuc ~ 


1 g cm~® (Hosokawa et _alj|2012a)). We obtain a radius of 


nuclear burning, which is given as 


3M \i8 M\¥ 
mue = ~ 1.2 17 . 
m (son ) ( Mo ) 


A comparison with Eq. (2) shows that this radius always 
remains smaller than the radius of the star by at least a 
factor of 1000, and in fact increases more gradually with 
stellar mass. A given mass shell M will only reach nuclear 
densities on timescales much longer than the initial Kelvin- 
Helmholtz timescale tg, which increases as R~!. We can 
thus determine the time when nuclear densities are reached 
by equating (QI) with the second term in (10). As a result, 
we obtain 


(11) 


(12) 


M 2/3 


During the time interval At, the star will accrete an addi- 
tional mass AM = MAt. Considering a given mass shell, 
we can calculate the ratio 

AM MAt ( M B 


a = A 710m 


13 
M M (1a) 
which describes how much additional mass is accreted be- 
fore the shell reaches nuclear densities. The ratio drops 
below 1 only for protostellar masses of 
M > 3.6 x 108 m3 Mo. (14) 
While for typical accretion rates of m ~ 1073, this hap- 
pens already at mass scales of order unity, this transition 
occurs only at substantially larger masses above 1000 Mo 
for rn > 1072. While a given shell evolves towards nuclear 
densities, the amount of additional mass that is accreted is 
thus substantially higher than the mass in that shell. From 
Eq. (12), we can further derive the mass in the nuclear core 
in the limit that t >> tinı( M). We obtain 


t 
Mnuc aad 
(= yr 


3/2 
) Mo = tro Mo, (15) 
where we defined t710 = t/(710 yr). As we neglected the 
term M/M in Eq. (22), we note that the derivation here 
implicitly assumes a sufficient mass supply to the protostar 
in order to feed the core. Clearly, this approximation will 
break down at the mass scale derived in (14). At that point, 
both the protostellar mass as well as the mass in the nuclear 
core may become approximately constant, and we expect a 
transition towards supermassive main sequence stars. 
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Fig. 2. The ratio of nuclear luminosity provided via the 


CNO-cycle over the Eddington luminosity of the star as a 
function of protostellar mass for different mass accretion rates. 
The CNO luminosity is calculated here as the minimum of the 
estimate in Eq. (24Jand the Eddington luminosity of the nuclear 
core (25). At high protostellar masses, one can clearly recognize 
the transition to the regime where nuclear burning is limited by 
the Eddington luminosity of the core. Larger accretion times 
imply that less time was available for contraction, thus reduced 
nuclear luminosities. 


From this model, we already obtain a set of relevant con- 
clusions concerning the evolution of the protostar. In par- 
ticular, if the accretion rate is high, a substantial amount 
of additional matter is accreted before a given shell reaches 
nuclear densities. As a result, the evolution of the protostar 
is dominated not by the interior shells, but by the additional 
matter which is accreted during the contraction. This be- 
havior will only change when the timescale of accretion be- 
comes comparable to the timescale on which the outer mass 
shell is able to reach nuclear densities. As we show here, 
the latter is significantly larger than the Kelvin-Helmholtz 
timescale of the star, and scales with the R7).. The tran- 
sition occurs when a critical mass scale of 3.6 x 10°73 Mo 
is reached, implying a transition towards a supermassive 
main sequence star. We note that the transition may even 
occur earlier if rn decreases with time. 


3. The nuclear evolution 


In this section, we will assess the potential role of nuclear 
burning. For this purpose, we show that the transition to 
the CNO-cycle rapidly occurs and regulates the production 
of helium in the nuclear core. We will then calculate un- 
der which conditions the fuel in the nuclear core will be 
exhausted while the protostar maintains its extended enve- 
lope. The latter implies the potential formation of a black 
hole in the interior, thus a transition into a quasi-star. 


3.1. Importance of the CNO-cycle 


While the pp-cycle will dominate in the very beginning (see 
appendix), we show here that the transition to the CNO- 
cycle rapidly occurs as a result of helium burning via the 


triple-a process. Following |Padmanabha (2000), the en- 
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ergy production rate via triple-a is given as 


5.1 x 1082Y ° 
in 7 P e7 4-4027/To erg gt s71, (16) 
9 


where Y ~ 0.25 denotes the mass fraction of helium and Tg 
denotes the gas temperature in units of 10° K. It is straight- 
forward to check that the energy production rate changes by 
several orders of magnitude for temperature changes of the 
order 10%. As shown by (2012a)), the high- 
est temperatures with Ty ~ 0.15 are expected in the center 
of the nuclear core, which has an enhanced nuclear density 
of 10 g cm~?. With a volume filling factor egu ~ 0.01 — 0.1, 
the contribution of that region to helium burning is sig- 
nificant due to the steep temperature dependence. For the 
following estimates, we will assume that efficient mixing oc- 
curs throughout the nuclear core, implying that the heavy 
elements produced in the central region will be available 
throughout the core. On the other hand, if mixing is in- 
efficient, nuclear burning via the CNO cycle will dominate 
even more in the central region due to its higher metallicity. 

Noting that the energy released by a single triple-a@ re- 
action corresponds to 1.166 x 10~° erg, the heavy element 
production rate per unit volume is given as 


nza = a Er. (17) 
The CNO mass in the core thus evolves as 

Meno = hsacnn Mnue X 12mp, (18) 
and an integration yields 

Mono = 7.1 x 107" epn t2 Mo. (19) 


Assuming efficient mixing, the resulting metallicity in the 
nuclear core is then 
_ Mono 


L= Maas = 6.3 x i07 t710 Efill- 


We recall that the energy production rate in the CNO cycle 
is given as (Padmanabha. 


Ta? 


(20) 


1, =i 


ergg. s (21) 


ECNO = 


Since the energy production rate depends sensitively on the 
temperature, we expect its contribution in the innermost 
core to be dominant. For the temperature Ty = 0.15, it 
is straightforward to show that energy production via the 
CNO cycle becomes comparable to the pp cycle for a metal- 
licity of Ze = 2 x 1071?. A comparison with Eq. (20) shows 
that the CNO cycle thus dominates after a short time of 
teno = 2.25€54 1 Y5, (22) 
where we introduced efu = 0.legu,—1. We note here that 
this timescale is probably not accurate, as our assumptions 
in section P] (in particular concerning the Eddington lu- 
minosity) only become valid at later times. Nevertheless, 
this result implies that the CNO cycle can be expected to 
be relevant early on, and thus needs to be considered. As 
the production of elements becomes increasingly efficient at 
late times, we expect that the expression (20) will yield a 
reasonable estimate in the period of interest. 


With the above assumptions, the CNO luminosity due 
to nuclear burning is given as 


Z 
Lono = €ono,Z. | = | Mnuc Est, (23) 


Ze 


where ecno,z. denotes the energy production rate via the 
CNO cycle evaluated at the critical metallicity Ze. Insert- 
ing (20), we obtain 


Loxo = 1.3 x 10°e2, #42 Lo. (24) 


The Eddington luminosity of the nuclear core is given as 


LFdd,core = 3.8 x 104112 Lo. 


(25) 
The ratio of these luminosities to the Eddington luminosity 
of the star is given in Fig. (2). Even here, the luminosity 
produced by nuclear burning approaches the Eddington lu- 
minosity of the star only around stellar masses of 10° Mo, 
implying no relevant impact on the stellar evolution during 
the earlier stages. A comparison of these expressions yields 
the timescale 


te = 2.1 x 10*eqt _yyr. (26) 


For e€gn,—1 ~ 1, as indicated by [Hosokawa et al] for 
m = 0.1, the CNO luminosity would exceed the Edding- 
ton luminosity of the core after ~ 2 x 10* yr, shortly after 
the end of their simulations. However, when the Edding- 
ton luminosity is reached, one would expect an expansion 
of the nuclear core, implying lower densities and an adia- 
batically decreasing temperature. Due to this thermostat, 
it is likely that the core will adjust to a state maintaining 
its Eddington luminosity. 

We will now demonstrate that the nuclear core is not 
converted to helium before the critical timescale te after 
which the core radiates at its Eddington luminosity. We 
consider the production of helium by the CNO process, 
given as 


Lono 


Nue,cNo = 43x 10° erg’ (27) 
and the helium mass production rate 

Mye,cno = 12my x Nue,cno- (28) 
An integration yields 

Mue,cno = 1.0 x 1073e2,, t7/2 Mo. (29) 


A comparison with the mass of the core indicates that they 
become comparable at 


tie,oNO = 2.25 x 10° eq 1 yr, (30) 


implying that a pure helium core will only be formed once 
the nuclear luminosity is equal to the Eddington luminosity 
of the core. 
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Fig. 3. The fraction of helium in the nuclear core, corre- 
sponding to the amount of exhausted fuel for nuclear burning, 
as a function of protostellar mass for different accretion rates. 
Higher accretion rates imply less time for protostellar contrac- 
tion. As a result, nuclear densities are reached later, implying 
less time for nuclear burning. 


3.2. Formation of a helium core 


We will now calculate the helium production rate under the 
assumption that the formation of a helium core occurs once 
the nuclear luminosity equals the Eddington luminosity of 
the core. Thus, 


Nue = iS (31) 
The helium mass production rate is thus given as 

Mye = 4mp x Nue (32) 
and an integration yields 

Mue = 1.05 x 1074 #34? Mo. (33) 


Under this assumption, the nuclear core becomes a pure he- 
lium core after ty ~ 6.8 x 10° yr. We note that at this point, 
the mass of the nuclear core is given as My ~ 9.4 x 10° Mo. 
In order to illustrate how the latter translate into total stel- 
lar masses, the ratio of helium mass to total mass in the 
core is displayed in Fig.B]for different accretion rates, show- 
ing that the hydrogen fuel is indeed exhausted for stellar 
masses of 3 x 10° — 3 x 10" Mo depending on the accretion 
rate. After that point, the burning of heavier elements will 
occur, but only last for a short period. As the mass of the 
core is considerably larger than the Chandrasekhar mass 
scale and the Tolman-Oppenheimer-Volkoff mass scale, the 
core may collapse and form a black hole within the star, a 
so-called quasi-star, as proposed by [Begelman et al] 
and (2010). The interior black hole may then 
continue accreting from the stellar envelope, while the ac- 
cretion rate of the quasi-star will only be limited by the 
Eddington rate corresponding to the total mass of the con- 
figuration. The derivation of tr and My however implicitly 
assumes that enough mass has been accreted in order to 
build up the core, requiring that 
My K tM. (34) 


Article number, page 6 of [9] 


In order to form a quasi-star, we thus require an accretion 
rate of M > My/ty, ie. 
m > 0.14. (35) 
For lower accretion rates, we expect that the extended en- 
velope of the star cannot be maintained at late times, im- 
plying that the star will contract and form a supermas- 
sive main sequence star. During this period, both the 
stellar mass and the mass in the nuclear core will be ap- 
proximately constant, implying a characteristic hydrogen 
burning timescale of 2 x 10° yr. For stellar masses above 
300 Ma, these are expected to collapse to massive black 


holes (Fryer et al] 2011). 


4. Discussion and outlook 


To investigate the interplay of accretion and Kelvin- 
Helmholtz-contraction of the protostar, we have presented 
an analytical model following the evolution of the mass 
shells from the extended envelope down to nuclear densi- 
ties. Based on this model, we are able to explain the funda- 
mental result by [Hosokawa et al] that the accreting 
protostar keeps its extended envelope beyond the adiabatic 
accretion phase. The reason is that the timescale on which a 
given mass shell may reach nuclear densities is much longer 
than the Kelvin-Helmholtz timescale of the star, as the pro- 
tostellar radius is considerably larger than the radius of the 
nuclear core. Instead of the Kelvin-Helmholtz timescale 
evaluated at the protostellar radius, the relevant timescale 
is thus the Kelvin-Helmholtz timescale at the expected core 
radius after contraction, providing an appropriate estimate 
of the evolutionary timescale. As a consequence, rapidly 
accreting protostars are expected to keep their extended 
envelopes until stellar masses of 3.6 x 108m Mo. At 
this point, the accretion timescale becomes longer than the 
Kelvin-Helmholtz timescale in the stellar interior, implying 
that the star will subsequently contract towards the main 
sequence. Then, we expect it to follow the typical evolution 
of a supermassive (proto)star (Shapiro & Teukolsky|[1986). 

We further assessed the impact of nuclear burning for 
the evolution of these stars. While deuterium burning will 
never be relevant during the protostellar evolution, hydro- 
gen burning via the pp-cycle could potentially affect the 
stellar evolution after 1015 years, which is however practi- 
cally irrelevant. However, due to the efficient metal pro- 
duction via the triple-a process, a transition to the CNO 
cycle is expected early on. The energy production by the 
CNO cycle is considerably more efficient, implying that the 
core can radiate at its Eddington luminosity. As a result, 
a helium core forms after 7 x 10° yr. The expected mass 
of the core is then 10® Mo, which provides an upper limit 
on the initial mass of the interior black hole. Therefore, if 
the available gas reservoir is sufficient to maintain accre- 
tion rates m >> 0.14 until this point, the central core may 
collapse during the protostellar evolution phase, and the 
object becomes a quasi-star. Alternatively, the protostar 
will evolve towards the main sequence to become a super- 
massive star. Such an object is then expected to collapse 


into a supermassive black hole at the end of its lifetime 
(Fryer & Heger|2011)). 


Our results have been derived based on a comparison of 
the mass accumulated by accretion with the characteristic 
mass scale where we expect a pure helium core. This mass 
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scale is based on the reasonable assumption that the lumi- 
nosity of the core is given by its Eddington luminosity, in 
agreement with the numerical results by 
(20124). Uncertainties are however present, as our cal- 
culation for the evolution of the core mass assumes that 
the luminosity of each mass shell is given by the Edding- 
‘ad mainte olka Coa. the 
latter is a good but not precise approximation. Further 
corrections can be expected in particular when approach- 
ing nuclear densities, when feedback from nuclear burn- 
ing becomes significant. Additional effects could be intro- 
duced as a result of rotation, which was not included in 
this model. We thus expect our results to provide an or- 
der of magnitude estimate for the critical accretion rate, 
which may be determined more accurately employing stel- 
lar evolution calculations, but also a more realistic model 
for the time-dependent accretion rate. We also note here 
that reported a slight dependence 
of their results on the employed boundary conditions for 
the protostar. While their main calculations adopted shock 
boundary conditions, they also explored the effect of lower- 
entropy accretion provided by photospheric boundary con- 
ditions. In that case, a slightly higher accretion rate of 
0.3 Mo yr! appears to be required to maintain the ex- 
tended envelopes, but the overall evolution remains very 
similar. 

For comparison, we note that recent simulations by 
Latif et_al! (2013d) reported accretion rates of ~ 1 Mọ yr7! 
in halos with ~ 10° Mọ. Adopting a constant accretion 
rate, the transition to a supermassive main sequence star 
would occur at a stellar mass of 3.6 x 10° Mo. As the nu- 
clear fuel is however exhausted earlier, a central black hole 
may form before this transition, giving rise to an evolu- 
tion as sketched by [Ball et al] (oii. However, due to the 
limited gas reservoir, and as here the halo mass is in fact 
comparable to the mass scale of the corresponding quasi- 
star, it seems likely that the accretion rate will substantially 
decrease at late times, implying that the transition towards 
the supermassive star should occur at an earlier stage which 
is then determined by the time evolution of the accretion 
rate, and that supermassive stars with 104 — 10° Me might 
be the most generic outcome of the collapse. More massive 
objects appear to be possible at least in principle, if a larger 
gas reservoir is available in more massive dark matter halos. 

In order to assess the potential influence of rotation, we 
estimated the amount of rotational support in the central 
103 Mo clumps reported by (20130, which var- 
ied between 5—20% in different simulations. With densities 
of ~ 1071? g cm~, these are still orders of magnitudes be- 
low the characteristic densities within the protostar, such 
that we cannot yet draw strong conclusions regarding the 
final amount of rotational energy. However, we note that 


reported a significant amount of rotation 
in primordial protostars based on numerical simulations of 
(2019), and a quite similar case can be expected 
here. As recently shown by(Reisswig et al] (2013), the latter 
may have a substantial impact on the collapse of supermas- 
sive stars, implying the potential formation of a black hole 
binary and a subsequent merger, accompanied with efficient 
emission of gravitational waves. The latter provides a po- 


tential pathway of probing black hole formation scenarios 
with LISA [] In a narrow mass range around ~ 55000 Mo, 


1 LISA webpage: |http://lisa.nasa.gov / 


one may further expect the occurence of highly energetic 
044 


supernovae with energies up to 10** erg, which can be po- 
tentially detected with JWSTE] (Johnsonet al]DOT). 
Returning to the fate of supermassive stars, their evo- 
lution in the presence of UV feedback was assessed by 
and in 
the case of spherical symmetry. These authors find that 
UV feedback is unable to stop accretion for rates above 
~ 0.1 Mo yr7!. As a result, the expected stellar mass is 


; ja \8/7 
Mov ~ 10°Mo | 5) - (36) 


For accretion rates of ~ 0.1 Mo yr™t, our results exceed this 


value, as the protostars remain on the Hayashi track up to a 
mass of ~ 3 x 10° Mo, making feedback inefficient. There- 
fore, taking protostellar evolution into account favours the 
formation of more massive objects. In all cases, a critical 
question concerns the time evolution of the accretion rate, 
since the transition to supermassive stars is regulated by 
the accretion rate at late times. While this paper provides 
a first assessment for the case of constant accretion rates 
and spherical symmetry, the potential implications of time- 
dependent accretion rates need to be addressed in more 
detail in the future, along with the implications of rotation 
during protostellar evolution. 
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Appendix A: Additional nuclear processes 


In section B.I] we have shown that the transition to the 
CNO cycle rapidly occurs, implying that the latter will reg- 
ulate the formation of a helium core as discussed in section 
[3.2] However, additional nuclear processes are expected to 
simultaneously occur, which we discuss here for complete- 
ness. In the two sub-sections below, we assess the role of 
deuterium shell burning as well as the impact of hydro- 
gen burning via the pp-change during the stellar evolution, 
showing that these processes will only have a minor impact 
on the stellar evolution. 


Appendix A.1: Deuterium burning 
As reported by [Hosokawa & Omukai) (2009), deuterium 


shell burning may occur even before nuclear burning starts 
in the core. The luminosity from deuterium shell burning 
is given as 


mh D/H 
Lp = 1.5 x 10°Lo (=) (5) 


where [D/H] denotes the deuterium abundance relative 
to hydrogen. For an accretion rate of m = 0.1, 
demonstrated that deuterium shell 
burning starts only at protostellar masses of ~ 80 Ma, and 
presumably even higher masses in case of higher accretion 
rates. A comparison with the Eddington accretion rate in 
Eq. (I) thus shows that the luminosity from deuterium shell 
burning will never become dominant in this mass range, but 
in fact becomes increasingly less relevant for larger masses. 

It is also clear that deuterium shell burning is not going 
to occur close to the outer surface. While the characteristic 
temperature for a given mass shell scales as T x GM/R, 


we note that, from Eq. (0), GM/R x t = M/M at 


times t > M/ M. In this regime, the central tempera- 
ture is almost spatially constant, consistent with the re- 
sults of (2012a), but already considerably 
larger than the deuterium burning temperature of ~ 10° K. 
Deuterium burning will thus occur when the first term in 
Eq. is still relevant. From Fig. [I it is evident that the 
radius of the shells changes considerably in that regime, 
while the mass in the shells is almost unchanged. Adopt- 
ing a scaling relation of T x GM/R with almost constant 
M implies that the radius has to change by 2-3 orders of 
magnitude for the atmospheric temperature of 5000 K to 
increase to a value above 10° K. We therefore expect that 
the radius of deuterium burning will correspond to a fixed 
fraction of the protostellar radius as long as the protostar 
maintains its bloated envelope. From a comparison with 
at a stellar mass of 1000 Mo, we 
obtain the normalization of this relation to be 

M 1/3 
Rp = 10 Ro (=) : 
As long as the protostar remains on the Hayashi track, the 
deuterium burning shell will not be able to move towards 
the atmosphere, and, therefore, it does not have an impact 
on the evolution of the protostar. 


(A.1) 


(A.2) 


Appendix A.2: The pp-cycle 


We consider first the energy production in the nuclear core 
via the pp-cycle, as the composition of the star is initially 


Dominik R.G. Schleicher et al.: Massive black hole factories 


0 
07 ji pae 
mally = 
#107 S 
A ee 
a 1076F eee =| 
= eo 
o7 H ee athe 4 
o genT m= 
oer m= 10. ------- 7 
0? pa iil sorora pop pil pa iil soror 1 Lini 
10? 10° 10* 10° 10° 107” 10° 
Protostellar mass [Moar] 
Fig. A.1. The ratio of nuclear luminosity provided via the 


pp-cycle over the Eddington luminosity of the star 
as a function of protostellar mass for different mass accretion 
rates. Larger accretion times imply that less time was available 
for contraction, thus reduced nuclear luminosities. 


primordial. As shown by|Hosokawa et_al (2012a), the typ- 
ical temperature in the core is ~ 108 K, and the nuclear 
3. We adopt here the expression of 


density ~ 1 g cm7 
Padmanabhar for the energy production rate, not- 
ing that 


2.4 x 104pX? _ 1/3 pe 
— e7338/T3" erg s7! go} 
T3 


(A.3) 


Epp = 


where Ty denotes the temperature in units of 10° K and 
X ~ 0.75 the mass fraction of hydrogen. While the pp- 
cycle considered here has only a moderate temperature 
dependence, we note that the CNO cycle scales as T?° 
around temperatures of 10° K, and the triple-a process as 
T*° around temperatures of 108 K. For the latter cases, 
we will thus need to take into account the higher tem- 
peratures within the center of the core, as they substan- 
tially contribute to the energy production via nuclear burn- 
ing. Now, with pnuc ~ 1 g cm~? and Ty = 0.1, we have 
€pp = 43 erg g`! s~'. The luminosity provided by the pp- 
cycle is then given as 


3/2 


Lpp = €ppMnuc = 21.5 5/6 Lo. (A.4) 


For comparison, the Eddington luminosity of the protostar 
is given as 


Lgaa = 3.8 x 104Mt Lo = 2.7 x 10" mtrioLo. (A.5) 


The ratio of these luminosities is given in Fig. (A.I), show- 
ing that it remains considerably smaller than unity for stel- 
lar masses up to 10° Mọ. Equating the two expressions, 
we find that the luminosity resulting from the pp-burning 
is relevant only at very late times tpp = 1.1 x 101° yr. 

For practical accretion rates of rn = 107° — 10 and pro- 
tostellar masses between 10° — 108 Mo, the contribution 
from the pp-chain will never be relevant, and can thus be 
neglected for the overall evolution of the stars. 

We now estimate the amount of helium produced in 
the core. During one fusion event, an average energy of 


4.3 x 107° erg is released, implying a helium production 
rate of 


Nye pp = m (A.6) 
The helium mass in the core thus evolves as 

Mite pp = 4Mp X Nee: (A.7) 
An integration yields 

Mue.pp = 6.0 x 1078 £2 Mo. (A.8) 


Equating the resulting mass with the total mass in the core, 
a helium core can be expected after a time of 1.2 x 10/0 yr. 
In summary, when only the pp-cycle is considered, nuclear 
burning is rather inefficient, and never going to have a sig- 
nificant impact on the evolution of the star. However, we 
will show in the next subsection that a transition to the 
CNO-cycle should be expected, and the implications of nu- 
clear burning will then become more relevant. 
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